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THEORY  OF  STATIC  FATIGUE  IN  BRITTLE  SOLIDS 
R.  H.  Doremus 
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Troy,  New  York  12181 


ABSTRACT 


A theory  of  static  fatigue  in  brittle  solids  is 
derived,  based  on  the  ideas  of  Hillig  and  Charles,  but 
using  a different  log  reaction  velocity- stress  relation- 
ship. The  resulting  equation  agrees  well  with  the 
functional  dependence  of  failure  time  on  stress  found 
experimentally,  and  experimental  values  of  the  stress- 
sensitivity  factor  are  consistent  with  the  theory. 
Variations  in  failure  times  with  sample  history  and 
glass  composition  can  only  partly  be  explained  by  the 
theory. 
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THEORY  OF  STATIC  FATIGUE  IK  BRITTLE  SOLIDS 


R.  H.  Doremus 

Rensselaer  Polytechnic  Institute 
Materials  Engineering  Department 
Troy,  New  York  12181 

INTRODUCTION 

Stress-accelerated  reaction  of  water  with  glassy  and  crystalline 
oxides  can  lead  to  weakening  and  crack  propagation.  This  phenomenon 
has  been  studied  in  two  different  ways.  The  time  to  fracture  of  a 
specimen  held  under  a constant  load  gives  a measure  of  "static 
iatigue",  and  the  rate  of  lengthening  of  macroscopic  (about  one  cm. 
long)  cracks  as  a function  of  load  demonstrates  slow  crack  propa- 
gation. in  some  materials  these  two  types  of  experiments  seem  to  be 
closely  related,  such  as  in  soda-lime  glass,  bui  in  other  materials 
such  as  fused  silica,  there  are  discrepancies  between  the  two  types 
of  measurements. 

Several  theories  of  static  fatigue  or  delayed  failure  of  brittle 

1 

solids  have  been  proposed;  that  of  Hillig  and  Charles  is  now 
generally  accepted.  Their  theory  considers  changes  in  the  dimensions 
of  cracks  ir,  a sample  surface  resulting  from  the  reaction  of  water 
with  the  material.  For  example,  in  a silicate  glass  the  load  is  born 
by  the  silicon-oxygen  network.  Water  breaks  silicon- oxygen  bonds  by 
the  reaction 

H20  + Si-O-Si  = SiOH  HOSi,  (1) 

so  this  reaction  can  change  the  dimensions  of  the  cracks  and  weaken 
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the  material. 

Defects  or  cracks  in  the  surface  of  a brittle  material  weaken  it 
because  stress  is  concentrated  at  the  crack  tip.  Consider  an  ellipti- 
cal crack  of  length  (semi-major  axis)  L and  half-width  (semi-minor 
axis)  a subjected  to  a tensile  stress  S perpendicular  to  L.  Let  x be 
the  distance  from  the  crack  axis  in  the  direction  of  the  applied 
stress;  x = 0 at  the  crack  tip  and  x = a at  the  half-width  of  the 

crack.  The  solution  to  elastic  deformation  of  the  crack  gives  a 

2 

stress  cr  at  the  crack  surface  distributed  as  follows  : 

g_  _ l+2L/a  - (L+a)  2 x2/a4  /2\ 

S 1+(l2  - a2)  x2/a4 

At  the  crack  tip  where  x = 0 

aQ  = S (1+2  L/a)  (3) 

j.f  L » a 

oQ  » 2S  L/a  = 2S  .Pl/Z  (4) 

where  p = a2/L  is  the  radius  of  curvature  of  the  crack  tip.  The 
stress  decreases  very  sharply  away  from  the  tip,  according  to  Eq.  2. 

Killig  and  Charles  assumed  that  the  rate  of  reaction  of  water 
with  the  material  controls  the  change  of  crack  dimension,  and  that 
this  rate  v is 

v = vc  exp  fa  (5) 

where  a is  the  tensile  stress  at  the  crack  tip,  j3  is  the  coefficient 
that  can  be  written  as  V*/RT,  where  V*  is  an  "activation  volume",  and 
vQ  is  the  rate  of  the  reaction  at  zero  applied  tensile  stress.  At 
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the  crack  tip  the  rate  at  which  the  crack  lengthens  dL/dt  = v,  so 
Eq.  4 also  gives  the  rate  of  slow  crack  propagation  as  a function  of 
tip  stress. 

In  order  to  relate  Eq.  (5)  to  the  time  of  delayed  failure,  the 
tip  stress  cr  must  be  expressed  as  a function  of  crack  dimensions 
(Eq.  4) . The  sharp  decrease  of  stress  away  from  the  crack  tip,  and 
the  exponential  relation  between  stress  and  reaction  rate,  mean  that 
the  tip  should  sharpen  much  more  rapidly  than  lengthen.  Thus  Hillig 
and  Charles  assumed  that  the  crack  length  L remains  nearly  constant, 
and  only  the  tip  radius  p changes  appreciably.  With  this  assumption 
they  solved  Eq.  5 in  terms  of  o.  The  crack  was  assumed  to  propagate 
to  failure  when  the  tip  stress  crD  equaled  ct,  the  ultimate  or  theo- 
retical cohesive  strength  of  the  material.  If  the  fracture  strength 
SN  at  liquid  nitrogen  temperature  is  time-independent,  reaction  (1) 
being  very  slow  at  this  temperature,  then  from  Eq.  4 

O't  = ^SN  f !,0/p o (6) 

where  Lo  and  pQ  are  the  initial  dimensions  before  corrosion.  The 
final  equation  for  fracture  time  t as  a function  of  stress  S wa.s 
found  to  be,  after  neglect  of  small  terms: 

In  t/tQ.5  = - Bcr  . (S/SN  - 0.5)  (7) 

where  tp  ^ is  the  fracture  time  at  S/SN  = 0.5.  The  linear  dependence 

of  log  fracture  time  on  stress  of  Eq.  7 is  found  experimentally  for 

3-5 

various  glasses  and  alumina  in  certain  stress  ranges  and  is  there- 
fore taken  as  support  of  the  theory. 
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There  is  some  question  about  the  assumption  that  L changes  much 
less  than  p.  Although  this  result  would  indeed  be  expected  from  Eqs. 
(2)  and  (5)  , as  soon  as  the  tip  has  corroded  a small  amount  the  stress 
distribution  at  the  tip  probably  is  not  given  by  Eq.  (2) . A tip 

O 

radius  of  about  20A  for  cracks  in  soda-lime  glass  was  found  in  an 
6 

etching  study  . If  the  applied  stress  is  0.25  SN,  then  the  tip  radius 

must  decrease  by  a factor  of  sixteen  for  failure  to  occur,  from  Eq.  4. 

This  gives  a final  tip  radius  of  less  than  atomic  dimensions,  which 

seems  unlikely  <,  Energy  considerations  also  indicate  that  this  small 

7 

a tip  radius  is  unlikely  . 

To  overcome  these  difficulties  a somewhat  different  approach  was 
8 

taken  . One  can  start  with  Eq.  1 in  the  form: 

v = °!li  = v exp  (1-a) 

dt  a 

The  assumption  is  made  that  a,  rather  than  I,,  is  constant  with  time, 
since  the  corrosion  rate  of  the  minor  axis  of  the  ellipse  is  much 
less  that  that  of  the  major  axis.  This  assumption  is  equivalent  to 
assuming  that  the  tip  radius  is  inversely  proportional  to  the  length 
during  corrosion,  since  p = a^/L.  With  this  assumption  Eq.  (7)  was 
derived,  predicting  a linear  dependence  of  log  t on  S/SN.  In  fact  a 
variety  of  assumptions  about  relative  changes  in  L and  o give  this 
dependence;  for  example,  if  it  is  assumed  that  p is  constant  and  only 
L changes,  a linear  relation  between  log  t and  S/SN  is  still  found. 
There  are  some  discrepancies  between  Eq.  7 and  experimental 
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results  on  static  fatigue  in  glasses.  The  factor  £crt  should  be 

invarient  for  a particular  glass,  since  it  includes  two  intrinsic 

parameters,  P , the  measure  of  the  effect  of  stress  on  a corrosion 

process,  and  the  ultimate  fracture  strength  of  the  material. 

Nevertheless  different  slopes  of  S/SN  versus  log  t plots  are  found 

for  different  surface  treatments  of  the  same  glass,  implying  difffer- 

ent  values  of  po  ^ for  these  treatments.  For  example,  in  etched 

soda-lime  glass  a value  of  77  is  found  for  /?ot,  compared  to  31  for 

3 

the  data  of  Mould  and  Southwick  . For  FN  borosilicate  glass  dif- 

4 

ferent  slopes  are  found  for  different  suriace  treatments  , and  for 

4 

fused  silica  fir  ^ is  31  for  abraded  samples  and  72  for  flame-p^.  j.xshed 
9 

samples  . 

Extensive  delayed  failure  tests  of  soda-lime  glass  after 

different  surface  treatments  and  an  analysis  of  distributions  of 

10 

failure  times  at  a particular  stress  showed  that  the  log  failure 

time  is  not  linear  with  applied  stress,  even  in  the  range  0.8  > 

S/SN  > 0.4.  Other  functional  dependencies  were  examined,  and  a good 

fit  to  the  data  was  found  for  inverse  proportionality  between  log 

failure  time  and  stress.  Such  a dependence  was  also  suggested  for 

experimental  delayed  failure  measurements  of  Baker  and  preston  on 
11,12  13 

Glass  , and  by  Taylor  „ Slow  crack  propagation  data  also  fit 

an  inverse  relation  between  log  crack  velocity  and  load  better  than 

10,14 

direct  proportionality  . In  this  paper  theoretical  reasons  for 
such  dependencies  are  explored,  and  the  equations  derived  are 
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compared  to  experimental  data  on  static  fatigue. 

DEPENDENCE  OF  REACTION  RATE  ON  STRESS 

The  effect  of  pressure  on  a chemical  reaction  rate  constant  k 

15 

is  found  from  transition  state  theory  to  be 


whfii’e  AV*  is  the  difference  between  the  molar  volume  of  the  "acti- 
vated complex"  of  the  reactants  and  their  normal  molar  volumes.  If 
AV*  is  constant  with  pressure, 

k = kQ  exp  (-  p AV*/RT)  . (9) 

In  the  equation  for  the  temperature  dependence  of  a chemical  reaction 
k = k1  exp  (-  Ah*/RT,  (10) 

where  AH*  is  the  activation  enthalpy  or  energy,  AH*  is  often  \nde~ 

pendent  of  temperature  over  wide  ranges  of  temperature.  However, 

16 , 17 

Av*  is  rarely  independent  of  pressure  so  Eq.  9 with  AV*  constant 
cannot  be  considered  a general  relation  for  pressure  or  stress 
effects.  Thus  the  assumption  of  Eq.  5 with  ? independent  of  stress 
has  no  experimental  support. 

An  equation  of  the  form 

v = exp  (-  cv/c ) (11) 

with  vro  and  ? as  stress-independent  coefficients,  is  proposed  for 
the  rate  of  reaction  of  water  with  oxides  subjected  to  a tensile 
stress  a;  is  the  limiting  reaction  rate  at  high  stress.  Equation 
11  with  v as  the  velocity  of  crack  propagation  and  cr  proportional  to 
the  load  fits  slow  crack  propagation  data  better  than  Eq.  5. 


L 


/ 
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Taylor  suggested  that  delayed  fracture  in  brittle  materials 
resulted  because  of  a time-dependent  rearrangement  of  the  structures 
of  the  material.  He  found  that  the  "activation  energy  for  fracture" 
should  be  reduced  as  the  bond  length  increased,  and  he  assumed  that 
water  also  reduced  the  activation  energy.  He  derived  an  equation  of 
the  same  form  as  Eq.  11  for  the  "rate  of  fracture".  These  ideas  can 
be  modified  somewhat  and  applied  to  the  present  point  of  viev;  to 
gain  some  possible  insight  into  the  origin  of  Eq.  11  and  parameters 
affecting  O'.  If  the  activation  energy  for  reaction  1 is  inversely 
proportional  to  the  change  in  bond  length  ha  from  the  original 
value  a0,  then 

£ = (12) 

0 RT£a  RTC 

where  B is  a constant,  so  o'  should  be  proportional  to  young's 
mo  lulus  E and  inversely  proportional  to  temperature.  The  constant  B 
should  be  proportional  to  the  activation  energy  for  the  reaction  of 
water  with  the  oxide  network. 


STATIC  FATIGUE 

If  the  reaction  of  water  with  an  oxide  as  a function  of  stress 
follows  Eq.  11,  then  from  Eq.  4: 


v 


DL 

dt 


= exp 


fc  a \ 
V2SL  / 


(13) 


Again  it  is  necessary  to  make  some  assumption  about  the  relative 

change  in  crack  length  L and  tip  radius  p as  corrosion  proceeds.  As 

2 

before  the  assumption  is  made  that  a is  constant,  or  that  p = a /L 
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is  inversely  proportional  to  L.  The  particular  relation  assumed 
does  not  change  the  final  equation  significantly,  as  shown  below. 
In  order  to  solve  Eq.  13,  the  following  substitution  is  made: 
&a 


X = 


2SL 


(-14) 


When  t = 0,  x = xQ  = 0:a/2SLo.  Then  Eq.  13  becomes: 

^ 4,=  - iS-Ss-i* 

x-^  cva 

or  integrating  between  t = 0 and  t = t: 


(15) 


efdz  _ 

2S 

V co  t 

x 

xo 

jral  is 

rx  ezd2  . 

[-ai 

X 

+ 

X 

1 1 

xo 

xo  x. 

„X  _Z 
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The  experimental  values  of  x are  several  times  larger  than  unity. 

Under  this  condition  the  last  two  integrals  can  be  expanded  as 
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follows  : 


px  j-f 


dz  = ( 1 + .-Li.  + iLi.  + ^4.  + 

■x-2 


•••) 


(18) 


z X ' X 

The  series  is  "semi-convergent",  that  is,  it  can  be  used  up  to  those 
terms  that  begin  to  increase  in  value.  If  x » 1,  then  the  approxi- 
mate result  for  Eq.  16  is: 

2S  v»  t 


exo 


f-a 


(19) 
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Since  xQ  = Of  S^/e^S,  when  S/Sjj  is  less  than  about  0.8,  ex  « ex°, 
and  finally 


exp 


(<*  SN  ) = ^ 


SN  t 


a 4-S 


'o  ^ t 


or  In  t 


= */gt 
S/S 


In 


N 


( SLofrt  \ 


(20) 

(21) 


Since  In  S/SN  varies  slowly  with  S/fn  for  S/SN  > 0.25,  this  equation 
gives  inverse  proportionality  between  log  of  the  fracture  time  and 
S/SN,  as  found  experimentally. 

Other  assumptions  about  the  relation  between  L and  p as  reaction 
proceeds  give  the  same  inverse  proportionality  between  In  t and  S/SN 
as  does  Eq.  21.  For  example,  if  it  is  assumed  that  c changes  much 
more  than  L,  as  suggested  by  Hillig  and  Charles,  the  exponent  of  x 
in  Eq.  14  approaches  one  instead  of  two,  and  the  final  equation  is 


In  t = + 

s/sN 


In  Lo 


(22) 


Alternatively  if  o is  assumed  to  be  constant  and  only  L changes,  the 
exponent  of  x in  Eq.  15  becomes  three,  and  with  x large  In  t is  still 
inversely  proportional  to  S/Sjj.  This  same  independence  of  the 
functional  relation  between  In  t and  S/SN  on  the  variations  of  L and 
p was  found  in  the  last  section  for  the  theory  of  Hillig  and  Charles. 
It  means  that  the  functional  dependence  of  log  t on  S/SN  (or  of 
crack  velocity  on  load)  ca  not  be  used  to  decide  between  different 
geometric  changes  in  the  crack  tip  and  length.  The  time  th  to 
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failure  at  S/SN  = h is,  from  Eq,  20: 

fch  = e*P  (?-  ) (23) 

STRUCTURE  OF  OXIDE  SURFACES 

In  this  section  recent  understanding  of  the  structure  of  oxide 
surfaces,  particularly  of  silicate  glasses,  is  summarized,  because 
this  structure  can  influence  the  rate  of  reaction  of  water  with 
the  oxidi 

A fracture  surface  of  an  oxide  contains  "dangling"  oxide  bonds. 
These  bonds  react  rapidly  with  atmospheric  water  to  form  metal 
hydroxyl  groups  such  as  SiOII,  BOH  or  A10II.  In  addition  to  these 
"external"  hydroxyl  groups  a hydrated  layer  several  atomic  distances 
thick  can  be  formed  on  the  oxide  surface  by  a reaction  like  that  in 
Eq.  1.  This  hydrated  layer  is  x>re  or  less  porous  to  water,  and  the 
interface  between  it  and  the  unhydrated  oxide  is  probably  the  actual 
load-bearing  point  Thus  a thin  crack  can  be  partially  filled  with 
a hydrated  oxide. 

The  surface  structure  of  a silicate  glass  containing  alkali  is 

even  more  complex.  Ion  exchange  between  alkali  ions  in  the  glass 

and  hydrogen  or  hydronium  ions  from  water  can  take  place : 

Na+  + 2H20  = H30+  + NaOH  (24) 

The  resulting  surface  structure  can  be  divided  into  layers,  as  shown 

19-21 

in  Fig.  1,  based  on  the  work  of  several  different  investigators 
If  the  glass  is  in  contact  with  liquid  water  the  outer  surface 
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dissolves  into  the  water;  in  contact  with  water  vapor  a layer  of 
alkali  hydroxide  or  carbonate  builds  up  on  the  glass  surface,  speeding 
corrosion  of  the  gel.  Movement  of  the  gel-glass  interface  results 
from  reaction  1,  and  the  rate  of  this  movement  can  be  associated  with 
the  velocity  v of  the  previous  section.  Since  the  mobility  of 
hyaronium  ions  is  much  smaller  than  that  of  sodium  or  lithium  ions, 
the  layer  containing  hydronium  ions  has  a much  higher  resistance  than 
the  unexchanged  glass.  On  the  other  hand,  the  gel  layer  is  quite 
open,  so  ions  in  it  have  a high  mobility.  These  differences  of 

21 

conductivity  of  the  various  layers  have  been  confirmed  by  Wilcby  in 
measurements  of  the  resistance  of  glass  membranes  progressively 
etched  by  HF. 

COMPARISON  OF  EQUATIONS  WITH  EXPERIMENTS 

In  Table  1 are  values  of  a/(?t  for  various  brittle  materials  as 

calculated  from  log  failure  time  versus  reciprocal  stress  plots  and 

Eq.  21.  All  these  plots  were  linear  for  the  data  in  Table  1,  are 

required  by  Eq.  21.  Soda-lime  samples  treated  in  a wide  variety  of 
3,10  4 

ways  all  came  close  to  the  same  rv/at  values.  Abraded  and 
9 

pristine  fused  silica  showed  widely  different  slopes  on  log  failure 

time  versus  stress  plots,  but  on  reciprocal  stress  plots  a*,  ut  the 

10 

same  slopes  and  C'/a-j.  values  were  found  . For  FN  borosilicate  glass 
the  values  of  o'/(Tt  was  dependent  on  treatment;  the  reason  for  this 
result  is  not  known. 

The  theoretical  fracture  stress  Cj.  can  be  calculated  from  the 


22 

following  equation  : 


12 


where  UQ  is  the  potential  energy  per  unit  surface  area  at  equili- 
brium separation  a,  and  E is  Young's  modulus.  UQ  can  be  found  from 
bond  energies  and  the  atomic  density  of  the  material.  A of 

1*3(10)^  Newtons/meter2  (2.6  (lO)^psi)  was  calculated  for  fused 
22 

silica  , which  is  not  much  greater  than  the  maximum  measured 
strength  of  fused  silica  at  -196°C,  or  1.4  N/m2  (2(l0)6psi).  For 
silicate  glasses  the  various  factors  in  Eq.  24  change  little  with 
glass  composition,  so  cr^.  should  be  nearly  the  same  for  most  of 
these  glasses.  A value  of  1.4 (10) 10N/m2  (2(10)6psi)  is  estimated 
for  soda-lime  glass. 

Composition  effects  on  c*  are  less  certain.  From  Eq.  12  and 
the  discussion  following  it  <y  should  not  change  much  for  different 
silicate  glasses,  since  Young's  modulus  for  most  of  these  glasses 
is  about  the  same,  and  reaction  1 governs  the  destruction  of  the 

23 

load-bearing  silicate  network  in  these  glasses.  Rana  and  Douglas 
measured  the  activation  energies  for  the  linear  corrosion  rates  of 
various  glass,  which  are  probably  related  to  the  activation  energies 
for  reaction  1,  as  described  in  the  last  section.  They  found  acti- 
vation energies  from  15  to  25  kcal/mole  for  various  soda-lime 
silicate  glasses,  so  it  might  be  reasonable  to  assume  a value  of 
20+5  kcal/mole . 
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Therefore  neither  & nor  should  change  much  for  silicate 
glasses,  in  agreement  with  the  similar  experimental  values  for 
Ot/crz  for  different  glasses. 

The  theoretical  fracture  strength  of  alumina  can  be  calculated 
from  Eq.  24  to  be  about  2.2 (10) 10N/m2  (3 . 2 (10) 6psi) , or  slightly 
larger  than  that  of  fused  silica.  The  activation  energy  for  the 
rate  of  reaction  of  water  with  aluminum- oxygen  bonds  is  uncertain, 
but  should  not  be  too  different  from  the  20  kcal/mole  for  silicates. 
Thus,  the  three-fold  increase  in  a/crt  found  experimentally  for 
alumina  as  compared  to  the  silicate  glasses  probably  results  mainly 
from  the  higher  modulus  of  alumina,  which  is  about  six  times  greater 
than  those  of  silicate  glasses. 

These  comparisons  show  that  the  values  of  <v/cr j-  found  experi- 
mentally are  in  reasonable  accord  with  theoretical  values,  and  lend 
support  to  the  present  derivation. 

Values  of  tn  for  a particular  material  treated  in  different 

ways  vary  widely.  For  soda-lime  glass  values  of  tn  differ  by  three 

3,10 

orders  of  magnitude  depending  upon  the  history  and  type  of  abrasion. 

The  tn  values  for  abraded  and  pristine  fused  silica  diverge  by  about 

4,9 

eight  orders  of  magnitude  . Between  alumina  samples  as-received 
and  samples  heated  through  several  cycles  up  to  about  1200 °C,  there 
was  a difference  of  about  two  orders  of  magnitude  in  tn.  These 
variations  cannot  be  readily  explained  by  Eq.  23. 

One  would  expect  «/<rt  for  a particular  material  to  be 


independent  of  treatment,  as  is  generally  found  experimentally, 
with  the  exceptions  noted  in  Table  1. 

Differences  in  initial  crack  length  LQ  can  be  deduced  from 
sample  strengths.  The  different  treatments  of  soda-lime  glasses 
led  to  variations  of  about  a factor  of  two  in  liquid  nitrogen 
strength.  Thus  the  initial  crack  lengths  should  vary  at  the  very 
most  by  an  order  of  magnitude,  since  the  radius  p of  the  crack  tip 
does  not  usually  change  much.  The  pristine  silica  had  a strength 
at  liquid  nitrogen  temperature  about  100  times  greater  than  the 
strength  of  abraded  fused  silica  at  this  temperature.  Thus  L0  for 
abraded  silica  could  be  four  or  five  orders  of  magnitude  greater 
than  for  pristine,  but  not  eight  orders  of  magnitude.  The  liquid 
nitrogen  strength  of  heated  alumina  was  about  1.5  times  that  for 
as- received  samples;  again  the  initial  crack  length  LQ  should  vary 
much  less  than  found  experimentally  for  tn. 

The  influence  of  various  treatments  on  the  reaction  rate  v0 
at  high  stress  is  uncertain.  One  would  expect  v^  to  be  propor- 
tional to  water  concentration  (humidity) . The  results  discussed 
here  were  carried  out  either  in  water  or  100%,  or  50%  relative 
humidity,  which  should  result  in  changes  of  no  more  than  a few 
times  in  vm.  For  fused  silica  and  alumina  v®  could  be  influenced 
by  the  thickness  of  hydrated  layer  on  the  oxide  surface,  but  the 
extent  of  such  influence  is  uncertain. 

The  layer  structure  of  surfaces  of  alkali  silicate  glasses, 
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described  in  the  last  section,  could  influence  vro  more  strongly. 

In  a very  thin  crack  near  its  tip  products  of  alkali-hydronium  ion 
exchange  (reaction  24)  could  build  up,  and  the  alkali  hydroxide 
increase  vro  substantially.  Various  treatments  could  influence  the 
basicity  at  the  tip,  changing  vm.  The  thickness  of  a pre-existing 
layer  of  hydrated  glass  in  the  crack  could  also  influence  the  rate 
of  ion  exchange,  thus  changing  the  alkalinity  and  vCJ.  The  extent 
of  such  influence  and  whether  or  not  it  can  explain  the  large 
experimental  changes  in  tn  await  more  detailed  experiments  on  the 
relation  between  the  structure  of  the  hydrated  layer  and  its 
dissolution  rate  in  water. 

With  the  large  uncertain  changes  in  tn  for  a particular 
material  it  would  seem  audacious  to  try  to  explain  differences  in 
tn  between  different  materials.  However,  for  certain  materials 
treated  in  the  same  way  there  seems  to  be  some  relation  between  tn 
and  the  chemical  durability  of  the  material.  In  Table  II  are 
listed  the  tn  values  for  three  glasses  and  alumina  that  had  been 
centerless  ground  or  abraded.  The  strengths  of  the  glasses  were 
all  about  the  same,  and  they  have  about  the  same  cy/Oj.  values  (see 
Table  I)  . ’ihus  torn  Eq.23  any  large  differences  in  tn  for  the  glasses 
should  result  from  differences  in  v,,,. 

Measurements  of  the  rate  of  dissolution  of  some  glasses  have 
been  made  by  various  authors,  and  some  of  these  results  are  given 


in  Table  III. 
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Douglas  and  co-workers  found  that  alkali- silicate  glasses 

dissolve  in  two  stages.  During  the  first  the  rate  of  dissolution 

is  proportional  to  the  square  root  of  time,  whereas  during  the 

second  the  rate  is  directly  proportional  to  time.  The  rate  in  the 

first  stage  is  controlled  by  interdiffusion  of  alkali  and  hydronium 

ions,  and  the  second  is  a steady- state  condition  in  which  the  rate 

19.  7 

is  controlled  by  dissolution  of  the  silicate  network  of  the  glass  . 

Thus  the  rate  of  dissolution  during  the  second  stage  corresponds  to 

the  movement  of  the  gel-glass  interface  described  in  the  last 

section  (see  Fig.  1) . For  sodium  and  potassium  silicate  glasses 

containing  a variety  of  third  oxides,  including  alkaline  earths, 

aluminum,  zirconium,  lead,  and  zinc,  the  rate  of  dissolution  during 

the  second  stage  at  100 °C  varied  by  no  more  than  a factor  of 
26 

twenty  , and  most  of  these  values  were  lower  than  that  for  soda- 
lime  glass  given  in  Table  Hi; 

From  the  data  in  Table  III  one  can  conclude  surprisingly  that 
the  rate  of  dissolution  of  silica  in  water  is  about  the  same  as 
that  of  alkali  silicate  glasses.  Therefore  if  vtt  shows  the  same 
changes  with  glass  composition  as  these  dissolution  rates,  the 
differences  in  tn  values  between  fused  silica  and  soda-lime  glass 
cannot  be  explained  by  changes  in  dissolution  rate.  There  is  an 
unexplained  difference  of  four  to  seven  orders  of  magnitude  between 
tn  for  silica  and  for  soda-lime  glass. 

Some  theories  for  the  rate  of  reaction  of  water  at  crack  tips 
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in  glasses  have  been  proposed  , but  they  do  not  take  into 
account  the  hydrated  layer  or  the  details  of  surface  structure 
given  in  the  last  section  and  Fig.  1.  Furthermore  they  cannot 
explain  the  experimental  variations  of  tn  with  humidity,  heat- 
treatment,  type  of  abrasion,  and  glass  composition. 

CONCLUSIONS 

The  present  theory  xs  claimed  to  be  an  improvement  over  that 
of  Hillig  and  Charles  because  it  removes  some  of  the  inconsistencies 
described  in  the  introduction.  The  functional  dependence  of  log 
failure  time  on  stress  in  static  fatigue  experiments  fits  Eq.  21 
derived  here  better  than  previous  equations  and  the  parameter  <y/rt 
found  experimentally  is  consistent  with  theoretical  expectations 
in  most  cases.  However,  Eq.  11  is  probably  not  valid  at  low 
stresses,  so  the  present  treatment  is  expected  to  break  down  at 
some  low  stress  value,  apparently  not  yet  achieved  in  static 
fatigue  experiments.  The  theory  in  its  present  form  does  not 
demonstrate  a fatigue  limit  at  low  stress?  since  soda-lime  glass 
samples  aged  in  water  become  stronger,  there  must  be  a fatigue 
limit  at  some  low  stress  for  this  glass.  However,  previous 
suggestions  of  fatigue  limits  in  static  fatigue  or  crack  propa- 
gation experiments  disappear  when  the  data  are  plotted  according 
to  Eqs.  11  or  21. 

Variations  in  tn  with  sample  history  or  between  different 
materials  can  only  partly  be  explained  by  the  present  theory. 
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Previous  theories  also  failed  in  this  respect.  Thus  unknown  factors 
change  the  sensitivity  of  samples  to  stress  corrosion  by  water,  even 
though  the  assumption  that  static  fatigue  results  from  stress- 
accelerated  reaction  of  water  with  the  material  seems  to  be  valid. 
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TABLE  I 

O'/or t FOR  DIFFERENT  MATERIALS,  CALCULATED  FROM 
STATIC  FATIGUE  MEASUREMENTS 


Materials 

Treatment 

£v/Cft 

Refs. 

Soda-lime  glass 

Abraded 

6.5 

3 , 10 

(Kimble  R-6) 

Etched 

11 

5 

Borosilicate  Glass 

Abraded 

5.8 

4 

(G.E.  FN, 

Centerless 

4 

Corning  7050) 

Ground 

11.8 

Fused  Silica 

Centerless 

) 

4,9,10 

(Amersil) 

Ground,  Pristine 

S 7.4 

96%  Alumina 

As-received 

•v 

(Wesgo  A1500) 

Ground  and  Fired 

} 21 

4 

Fired  only 

V 

\ 


II 


TABLE  II 


LOG  FAILURE  TIME  IN  SECONDS  AT  S/SN  =r  0.5  FOR  VARIOUS 
MATERIALS  CENTERLESS  GROUND  OR  ABRADED 


Material 

Treatment 

Log  fcn 

Liquid  Nitrogen 
. Strength  SN 
MN/m2  (K  psi) 

Ref. 

Alumina 

Centerless  Ground 

4.8 

520  (72) 

4 

Fused  Silica 

II 

4.9 

100  (15) 

4 

Fi'T  Borosilicate 

► M 

1.5 

100  (15) 

4 

glass 
Soda- lime 

Various  abrasions, 

1.0  to 

70  to 

V.O 

!— * 
O 

Silicate  glass 

• aged  in  water 

-2.0 

150  (10-22) 

2,0 


tr< 


TABLE  III 

RATES  OF  DISSOLUTION  OF  VARIOUS  GLASSES  AT  25 °C 


Glass 

Medium 

Silica 

0.9%  Nad 
0.1%  NaHC03 

1M  NaCl04 
in  v;ater 

Lithium  Silicate 
«25%  Li20  plus 
other  oxides 

water 

Soda- lime  Silicate 
15%  Na20f  10%  CaO 

water 

Rate  of  Dissolution 


cm/sec . 

Ref 

4(10)“13 

24 

3 (10) ~12 

25 

io-13 

20 

2(lU)-12 

23 

(extrapolated) 
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Fig.  1 Schematic  structure  of  the  surface  of  an  alkali  silicate  glass. 
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CRACKS  AND  ENERGY  - CRITERIA  FOR  BRITTLE  FRACTURE 

R.  H.  Doremus 

Rensselaer  Polytechnic  Institute 
Materials  Engineering  Department 
Troy,  New  York  12181 

! INTRODUCTION 

j 

* . Brittle  solids  fail  by  propagation  of  flaws,  usually  in  their 

surfaces.  These  flaws  are  thought  of  as  thin,  deep  cracks  pene- 
trating from  the  surface  into  the  solid.  At  the  tip  of  such  cracks 

i 

any  applied  tensile  stress  becomes  magnified,  leading  to  propagation 

of  the  cracks  at  applied  stresses  S much  lower  than  the  ultimate  or 

theoretical  cohesive  stress  crt  needed  to  break  apart  the  solid. 

1 

Inglis  solved  the  elastic  equations  for  an  elliptical  crack 
subjected  to  a tensile  stress  S perpendicular  to  tte  crack  and 
found  the  stress  a at  the  crack  surface  to  be: 

(1) 

where  c is  the  depth  (semi-major  axis)  of  the  crack,  a is  the  half- 
width (semi-minor  axis) , and  X is  the  distance  from  the  crack  axis. 
At  the  crack  tip  X = 0,  and: 

0-o  = S(1  + 2c/a)  . (2) 

If  the  crack  is  deep  and  thin  c » a,  and: 

crQ  « 2S  c/a  = 2S  /c/o  (3) 

where  p = a2/c  is  the  radius  of  curvature  at  the  crack  tip.  Even 

1 

if  the  crack  is  not  elliptical  the  tip  stress  is  given  by  2S  /c/p. 
This  relation  shows  that  the  stress  at  the  tap  of  a deep,  thin 


21  = l+2c/a  - (c+a)2  X2/a^~ 
s l+(c2-a2)  x2'/a4 


2 


crack  is  much  greater  than  the  applied  tensile  stress;  this  magnifi- 
cation explains  the  low  fracture  stresses  of  brittle  solids  compared 
to  their  theoretical  cohesive  strengths  tf-t.  The  fracture  stress 
Sf  is 


Sf  = o ‘t  ^.c/  3/"c« 


(4) 


. Griffith  calculated  the  strain  energy  £ per  unit  length  of 
an  elliptical  crack  of  depth  c subjected  to  a tensile  stress  S and 
found  that  the  crack  lowered  the  energy  of  a stressed  specimen  by 


an  amount,  in  plane  stress: 


£ = - 


*c2S2 


(5) 


where  Y is  Young's  modulus  for  the  material.  Griffith  assumed  that 
the  elastic  energy  released  during  crack  propagation  provides  the 
additional  surface  energy  of  the  growing  crack,  4yL  per  unit  crack 
length,  where  y is  the  surface  energy,  and  derived  the  following 


equation  (as  modified  by  others)  by  equating  the  differentials  of 


the  two  energies  with  respect  to  length: 


si  = ZLV 


rr  c (6) 

where  Sf  is  the  applied  tens  .le  stress  needed  to  propagate  a crack 
or  depth  c. 

A number  of  subsequent  authors  have  used  energy  criteria  to 
understand  fracture  and  crack  propagation.  However,  Eqs.  4 and  6 


both  profess  to  serve  as  criteria  for  fracture,  and  yet  they  are 

inconsistent.  Various  ways  of  reconciling  these  inconsistencies 

3,4 

have  been  suggested 


3 


In  this  paper  energies  of  crack  propagation  are  considered  in 

terms  of  the  first  and  second  laws  of  thermodynamics,  and  it  is 

concluded  that  the  Griffith  criterion  provides  a necessary  but  not 

sufficient  condition  for  crack  propagation  and  failure,  as  suggested 
5 3 

by  Hillig  . An  argument  to  the  contrary  is  considered,  and  surface 
energies  of  brittle  solids,  particularly  oxides,  are  briefly 
discussed.  Implications  of  the  present  view  to  experiments  and 
theories  on  fracture  are  appraised. 

FIRST  LAW  OF  THERMODYNAMICS 

The  system  under  consideration  is  a block  of  material  in  a much 
larger  solid  subject  to  a uniform  tensile  stress  S.  One  surface  of 
the  solid  contains  a thin  crack  of  uniform  dimensions  along  its 
length,  which  passes  through  the  system  or  block.  The  crack  extends 
much  beyo.\d  the  ends  of  the  block,  so  there  are  no  end  effects.  The 
outer  dimensions  of  the  block  are  assumed  to  remain  constant  with 
small  changes  in  crack  dimensions.  The  length  of  the  crack  from  one 
face  to  the  other  is  one  length  unit,  so  all  energies  are  for  unit 
crack  length. 

The  first  lav/  of  thermodynamics  states  that  for  a thermodynamic 
system  there  exists  a state  function  E called  internal  energy  such 
that 

dE  = dq  - dw  (7) 

where  dq  is  the  heat  added  to  the  system  from  the  surroundings  and 
dw  is  the  work  done  on  the  surroundings  by  the  system. 

Temperature  changes  are  usually  negligible  in  elasticity 

27 


4 


6 

problems.  As  Timoshenko  and  Goodier  point  out  , the  difference 
between  the  adiabatic  and  isothermal  moduli  of  elasticity  are  small - 
for  iron  the  difference  is  0.26%.  Thus  temperature  changes  in  crack 
propagation  problems  are  often  ignored,  and  only  mechanical  and 
surface  energies  considered.  I contend  that  temperature  changes 
during  crack  propagation  cannot  be  neglected,  and  that  they  play  an 
important  part  in  the  energy  balance  of  Eq.  7.  There  is  no  direct 
experimental  evidence  to  support  or  refute  this  contention,  in  an 

attempt  to  measure  the  temperature  change  in  glass  near  a propagatir.. 

7 

crack  , a thermister  placed  a few  tenths  of  a millimeter  from  the 
crack  path  showed  no  change  in  temperature  to  within  an  accuracy  of 
a few  tenths  of  a degree.  However,  temperature  changes  considerably 
smaller  than  this  would  still  play  an  important  role  in  Eg.  7. 

Two  extreme  possibilities  are  that  crack  propagation  occurs 
isothermally  or  adiabatically.  If  the  process  is  to  be  isothermal, 
heat  must  be  added  to  or  removed  from  the  surroundings.  Since  the 
crack  propagates  rapidly  it  seems  more  realistic  to  consider  an 
adiabatic  process  where  no  heat  is  exchanged  with  the  surroundings 
(dq  = 0)  and  the  temperature  of  the  material  can  change.  In  either 
case  if  the  system  eventually  reaches  the  same  state  the  energy 
change  dE  will  be  the  same,  since  E is  a state  function. 

Consider  the  changes  in  Eq.  7 as  the  depth  c of  the  crack  chang. 
adiabatically  by  dc.  No  external  work  is  done  on  the  system  by  the 
surroundings,  because  the  dimensions  of  the  block  remain  cr nstant 

as 


5 


and  there  is  no  normal  stress  on  the  crack  surface.  Thus  both  dq  and 
dw  are  zero,  and  all  energy  changes  are  changes  in  dE.  The  stored 
elastic  energy  £ decreases  by  an  amount: 
d£  - 2tt  cS2  dc 

Y (8) 


from  Eq.  5.  The  surface  increases  by  an  amount  2dc  per  unit  crack 
length,  so  the  increase  in  surface  energy  is  4y  dc.  The  temperature 
of  the  material  near  the  crack  can  increase  or  decrease;  the  amount 
of  heat  involved  in  this  change  is  dQ,  not  to  be  confused  with  dq, 
the  heat  exchanged  with  the  surroundings.  Then  F~  n becomes: 


2^  cS2dc  , 

dQ  - + 47  dc  = 0 


(9) 


Only  if  temperature  changes  are  negligible  (dQ  = 0)  can  the  Griffith 
equation,  Eq.  6,  be  derived  from  Eq.  9.  If  the  fracture  stress  is 
larger  than  predicted  by  Eq.  6,  then  cS2  dc/Y  > 4y  dc,  dQ  is 
positive,  and  the  temperature  of  the  material  near  the  crack 
increases.  If  2ir  cS2  dc  <47  dc,  then  the  temperature  decreases. 

It  is  interesting  to  estimate  the  temperature  change  expected 
from  a propagating  crack.  Consider  a crack  for  which  the  Griffith 

O A 

fracture  stress  (Eq.  6)  would  be  70MN/m  (10  psi.) . If  the  actual 
fracture  stress  is  double  this  value,  then  about  20  joules  of  energy 

per  cm2  of  crack  surface  formed  are  available  for  heating  the  sample. 

•3  °C 

If  the  density  is  2.5gms/cnr  and  the  heat  capacity  800  joules/gm  , 

a region  half  a millimeter  thick  at  the  crack  surface  will  be  heated 


about  0.1°C. 


6 
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The  conclusion  from  this  discussion  is  that  a criterion  for 
fracture  can  be  derived  from  an  energy  balance  alone  only  by 
assuming  negligible  temperature  change,  for  which  assumption  there 
is  no  basis. 

SECOND  LAW  OF  THERMODYNAMICS 

The  second  lav;  of  thermodynamics  can  be  stated  for  the  present 
system  as  fallows:  For  a single  phase  there  exists  a state  function 
ri  called  entropy  such  that 

dE  = - PdV  + T dr?  + 7 dA  + d£  (10) 

when  no  mass  is  added  to  or  subtracted  from  the  system,  and  P is  the 
pressure  and  V the  volume  of  the  system.  A second  part  of  the 
statement  is  that  for  any  body  consisting  of  one  or  more  phases  and 
isolated  from  the  surroundings  (E  and  V constant) , the  entropy  of 
the  body  increases  or  remains  constant  according  to  whether  the 
body  is  or  is  not  in  internal  equilibrium: 

dr?  s 0 (E  and  V constant)  (11) 

For  the  present  system  V is  constant,  so  from  Eq.  10: 

T dr)  = dE  - IdA  - d&  (12) 

Then  f i om  Eq.  (9) 

T dt;  = dQ  (13) 

If  the  stress  is  lower  than  the  fracture  stress  of  Eq.  6,  then  the 
temperature  of  the  material  around  the  crack  decreases  and  dQ  is 
negative.  However,  this  condition  is  contiary  to  Eq.  11.  Thus 
crack  propagation  cannot  occur  for  stresses  below  that  predicted  by 
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Eq.  6,  and  this  equation  provides  a necessary  condition  for  crack 

propagation.  This  requirement  based  on  the  second  law  of  thermo- 

3 

dynamics  was  also  recognized  by  Orowan  . 

It  is  interesting  to  inquire  about  what  happens  physically 
when  the  applied  stress  is  high  enough  to  cause  propagation  by  Eq.4# 
but  propagation  cannot  occur  because  is  lower  than  calculated 
from  Eq.  6.  The  stress  at  the  crack  tip  will  be  higher  than  the 
theoretical  ultimate  strength  of  the  material;  the  material  at 
the  tip  will  fluctuate  to  a larger  tip  radius  so  that  the  tip  stress 
is  reduced  below  <rt. 

The  minimum  tip  radius  Pm  for  which  fracture  will  occur  can  be 
found  by  equating  Eqs.  4 and  6; 


Pm  " 


by  y 

^7 


(14) 


For  fused  silica,  a representative  brittle  solid,  Y is  7.2(10)  ^ 

N/m^,  and  theoretical  calculations  give  Crfc  = 1.8  (10)  ^M/m2  and 

7 =2.9  J/m^.  Then  from  Eq.  14  Pm  is  about  16A. 

This  discussion  shows  that  any  theory  of  fracture  in  which  the 

tip  radius  of  a crack  is  less  than  pm  cannot  be  correct.  In  parti- 

8,9 

cular  some  theories  have  postulated  a cusp- shaped  tip,  with  an 
"infinite"  radius  of  curvature.  As  soon  as  a tensile  stress  large 
enough  to  give  a tip  stress  greater  than  is  applied  to  such  a 
crack,  it  fluctuates  to  a larger  tip  radius  without  propagating. 
Thus  such  a tip  shape  cannot  be  significant  in  crack  propagation. 


t- 


8 


AN  ARGUMENT  TO  THE  CONTRARY 
3 

Orowan  argued  that  the  Griffith  criterion  of  Eq.  6 is  both 
necessary  and  sufficient,  and  tried  to  establish  necessity  from  an 
atomic  model.  The  present  arguments  are  entirely  macroscopic,  with 
no  recourse  to  an  atomic  model.  Orowan  assumed  an  intermolecular 
stress-distance  relation  illustrated  in  Fig.  1.  At  low  stress  the 
stress-displacement  relation  is  linear  (Hooke's  law);  as  the 
ultimate  fracture  stress  is  approached  the  stress  for  unit  dis- 
placement decreases.  Orowan  assumed  that  there  is  a region  near 
the  crack  tip  where  the  stress  is  constant  because  of  the  failure 
of  Hooke's  law  at  high  stress.  Neuber  showed  that  if  such  a region 
exists,  the  stress  in  it  will  be  about  the  same  as  at  the  tip  of  a 
crack  of  the  same  length  and  with  a tip  radius  about  the  same  as 
the  radius  of  the  constant-stress  region.  If  this  radius  is  about 
equal  to  the  atomic  spacing  b,  then  from  Eq.  (4)  the  fracture 
stress  is 

Sf  = crt  / b/2  / c (15) 

Various  equations  for  f*t  have  been  derived;  Orowan  used: 

Ct  = /2Y  y/b  (16) 

If  this  value  is  substituted  into  Eq.  (15)  the  result  is 

Sf  = /Yy/2c  (17) 

which  is  nearly  the  same  as  Eq.  (6) . Thus  Orowan  claimed  that  if 
the  applied  stress  is  the  value  given  by  Eq.  (6) , the  stress  at  the 
tip  of  the  crack  will  be  ?t,  and  failure  will  occur. 

32- 


Consider  a crack  of  depth  c with  a tip  radius  p » pm;  in 
other  words,  a crack  with  a tip  radius  much  greater  than  the  value 
from  Eg.  14  for  the  Griffith  criterion  to  aj.ply.  Now  increase  the 
tensile  stress  S on  this  crack.  When  S reaches  the  value  calculated 
from  Eq.  6 the  stress  at  the  tip  will  still  be  much  less  than 
so  Hooke's  law  will  hold  around  the  crack  tip.  Only  when  the 
applied  stress  reaches  a much  higher  value  than  calculated  from 
Eq.  6 do  deviations  from  Hooke's  law  appear.  Thus  the  propagation 
stress  for  this  crack  is  greater  than  predicted  by  Eq.  6,  and 
Griffith's  criterion  fails.  Orowan's  argument  applies  only  to 
cracks  with  tip  radii  smaller  than  pm. 

There  is  no  experimental  support  for  the  stress-displacement 

curve  shown  in  Fig.  1,  although  it  is  often  used  for  theoretical 

5 

calculations.  Hillig  measured  the  stress-strain  behavior  of  fused 
silica  up  to  stresses  close  to  the  ultimate  theoretical  stress;  the 
elastic  modulus  actually  increased  slightly  at  the  highest  stresses, 
and  there  was  no  indication  of  any  decrease  as  predicted  by  the 
curve  in  Fig.  1.  Thus  there  is  no  reason  to  expect  a constant- 
stress  region  at  the  crack  tip  as  postulated  by  Orowan. 

SURFACE  ENERGY 

The  surface  energy  y used  here  may  be  considered  to  be  one- 
half  the  energy  needed  to  separate  two  atomic  planes  of  unit  area; 
it  can  be  approximately  calculated  from  th  bond  energy  and  atomic 
density  of  a solid.  In  this  way  y values  of  2.9  joules/m2 
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2 2 
(2900  ergs/cm  ) for  fused  silica  and  1.2  joules/m  for  cv-alumina 

can  be  calculated. 

In  an  actual  fracture  experiment  surfaces  formed  may  react 
with  ambient  gases.  Freshly  fractured  oxide  surfaces  are  particu- 
larly reactive;  they  form  hydroxyl  groups  on  their  surfaces  by 
reaction  with  water,  for  example; 

Si-0-  + Si-  + H20  = 2Si0H  (18) 

The  surface  energy  of  an  oxide  surface  covered  with  hydroxyl  groups 
is  much  lower  than  calculated  above.  This  is  the  reason  that 
measured  values  of  surface  energies  of  molten  oxides  are  much 

lower  than  the  calculated  values. 

11 

Nadeau  fractured  glass  cylinders  at  high  stresses  in 

compression  into  fine  fragments,  and  measured  the  heat  given  off 

during  fracture  calorimetrically . The  amount  of  heat  given  off  was 

not  much  less  than  the  elastic  energy  stored  in  the  glass  just 

before  fracture.  It  seems  likely  that  this  heat  results  from  the 

exothermic  reaction  18  occurring  following  fracture. 

Surface  energies  have  been  calculated  from  fracture  measure- 
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ments  in  several  different  ways  . In  all  of  these  methods  the 
Griffith  relation  (Eq.  6)  is  assumed  valid.  In  one  method  the 
stress  needed  to  fracture  a material  with  a crack  of  depth  c is 

measured,  and  the  surface  energy  7 is  calculated  from  Eq.  6; 

c 

7 = — 


2Y 


(19) 


In  the  present  view  the  stress  needed  to  propagate  a crack  of  length 

c will  be  greater  than  or  equal  to  the  value  calculated  from  Eq.  6. 

Thus  one  might  expect  values  of  Y calculated  by  these  methods  to  be 

larger  than  the  theoretical  values  and  scattered.  These  are  the 

16 

results  found  experimentally.  Y values  for  sapphire  range  from 

14,16 

6.0  to  32  j/m  ,♦  for  polycrystalline  alumina  , from  18  to 
2 15,17 

54  j/m  ; for  soda-lime  glass  , from  3.4  to  11  j/m2;  for 
13 

magnesia  from  1.1  to  3.0  j/m2.  It  seems  reasonable  to  conclude 
that  the  lowest  value  for  each  material  gives  only  an  upper  bound 
to  the  "true"  surface  energy. 

DISCUSSION 

Brittle  solids  fail  because  of  the  magnification  of  stress  at 
the  tips  of  thin,  deep  cracks  (Eq.  3) . Thus  it  is  not  surprising 
that  Eq.  6,  in  which  this  magnification  plays  no  role,  cannot 
provide  a criterion  for  fracture.  In  the  integration  to  calculate 
the  change  in  elastic  energy  caused  by  the  presence  of  the  crack 
(Eq.  5)  the  terms  involving  the  high  stress  at  the  crack  tips  are 
negligible,  because  the  volume  over  which  this  high  tip  stress 
occurs  is  very  small,  and  becomes  smaller  as  the  tip  stress  increases 
In  a brittle  material  some  surface  cracks  can  have  the  minimum 
radius  pm  given  by  Fq.  14,  because  cracks  starting  with  smaller 
radii  will  fluctuate  to  this  value  as  soon  as  the  required  stress  is 
applied.  Thus  in  many  instances  brittle  materials  with  cracks  of 
depth  c will  fracture  at  the  minimum  scress  as  calculated  from  Eq.b. 


12 
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In  some  experiments  on  fracture  of  soda-lime  glass  the  crack  depth 

was  measured  from  etching  experiments  and  the  tip  radius  calculated 

to  be  about  20A.  This  is  close  to  the  radius  calculated  from  Eg. 

14,  so  the  Griffith  equation  (Eq.  6)  should  apply  in  this  case. 

However,  when  this  glass  was  treated  in  water  its  strength  increased, 

19 

yet  the  crack  length  remained  the  same  . In  this  case  the  Griffith 
equation  was  not  obeyed,  and  the  tip  radius  was  calculated  to  be 

O 

33A,  or  considerably  greater  than  the  value  from  Eq.  14.  These 
results  emph  size  that  in  fracture  experiments  the  tip  radius  of 
the  crack  is  a parameter  that  cannot  be  ignored,  as  is  done  in  an 
energy  approach;  Eq.  3 provides  the  correct  criterion  for  crack 
propagation  and  failure,  not  Eq.  6. 
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